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Abstract 

We present preliminary results of a search for several exclusive charmless hadronic B decays from 
electron-positron annihilation data collected by the BABAR detector near the T{4S) resonance. 
These include three-body decay modes with final states h^h^h^ and /i^/i^vr", and quasi-two-body 
decay modes with final states X^h and X^K^, where h = tt or K and X^ = rj' or u. We find 
B{B^ p^Fvr^t) = (49±13l^) x 10^^ and B{B+ r]'K+) = (62 ±18 ±8) x 10"^ and present upper 
limits for eight other decays. 
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1 Introduction 



The charmless hadronic decays of neutral B mesons of interest in this paper can be used to explore 
CP violation arising from two possible interference effects. Indirect CP violation in decays 
can arise from interference between the direct amplitude and one involving B^-B^ mixing, while 
direct CP violation in charged B decays can occur via interference between tree and penguin 
amplitudes. The former offers the future possibility of measuring directly the CKM angle a of the 
Standard Model. It is expected that this will require samples of order 100 million B mesons. 
The magnitude of direct CP violation is more difficult to estimate but interesting constraints can 
be provided with much smaller samples. In addition, the tree-penguin interference may allow 
measurement of the CKM angle 7 [§] , via measurements of the decay rates of modes presented here 
and the related B — > Kir and i? — > tttt decay modes Q. 

The initial challenge is to obtain significant samples of these rare decays. Many of the modes 
of interest have only recently been observed for the first time, or remain undiscovered [Q, |5|, |^. 
In this paper, we describe preliminary searches of the initial BABAR data sample for a number of 
charmless hadronic B decays, and give preliminary measurements of their branching fractions. The 
three-body final states are limited to those with at least two charged tracks (vr^ or K^) and at 
most one vr'^. The hadronic resonances p and K* are sufficiently short-lived that the modes which 
contain them interfere quantum mechanically with other relevant 3-body final state modes. We 
have however searched for them independently using appropriate kinematic selections within the 
Dalitz plot. The quasi-two-body decays involve an rj' or uj meson accompanied by a neutral kaon 
or charged pion or kaon. We summarize the decay modes considered^ as follows: 

three-body quasi-two-body 



B+ ^ p^K+ ^ 



s 



B+ K+'K~TT+ B + 

B+ p°7T+ B+ ujK+ 

B+ 7r+7r-7r+ B^ ujK^ 
B^ ^ p^vr^ 



2 The BaBAR detector and data 

The data used in the analyses were collected with the BABMi detector at the PEP-II storage 
ring. The BABMi detector, described in detail elsewhere 0, consists of five active sub-detectors. 
Surrounding the beam-pipe is a silicon vertex tracker (SVT) to track particles of momentum less 
than ~120MeV/c and to provide precision measurements of the positions of charged particles 
of all momenta as they leave the interaction point. A beam-support tube surrounds the SVT. 
Outside this is a 40-layer drift chamber (DCH), filled with an 80:20 helium-isobutane gas mixture 
to minimize multiple scattering, providing measurements of track momenta in a 1.5 T magnetic 
field. It also provides dB/dx measurements to help charged particle identification. Surrounding 
the outer circumference of the drift chamber is a novel detector of internally refiected Cherenkov 
radiation (DIRC) which provides charged hadron identification in the barrel region. This consists of 

^Charge conjugate decay modes are assumed throughout this paper. 
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quartz bars of refractive index ~1.42 in which Cherenkov hght is produced by relativistic charged 
particles. This is internahy reflected and collected by an array of photomultiplier tubes, which 
enable Cherenkov rings to be reconstructed and associated with the charged tracks in the DCH, 
providing a measurement of particle velocities. Outside the DIRC is a CsI(Tl) electromagnetic 
calorimeter (EMC) which is used to detect photons and neutral hadrons, and to provide electron 
identification. The EMC is surrounded by a superconducting coil which provides the magnetic field 
for tracking. Outside the coil, the flux return is instrumented with resistive plate chambers (IFR), 
interspersed with iron which may be used for the identification of muons and mesons. 

The data sample used for the analyses contains 8.8 million BB pairs Q, corresponding to 
7.7fb-i taken on the T{AS) resonance. In addition, 1.2 fb ^ of data taken off-resonance have been 
used to validate the contribution to backgrounds resulting from e^e~ annihilation into light qq 
pairs. These data have all been processed with reconstruction software to determine the three- 
momenta and positions of charged tracks and the energies and positions of photons and merged 
TT^ s. Refined information on particle type from the various sub-detectors described above is also 
provided, and is used in particle identification algorithms in the analyses, as described below. 

3 Candidate selection 

Charged tracks are required to satisfy some standard track criteria, including a requirement that 
the momentum is less than lOGeV/c and that the transverse momentum is greater than 0.1 GeV/c. 
They are required to have at least 20 hits in the DCH and to originate close to the beam-spot. 

Photon candidates are identified in our calorimeter as deposits of energy unassociated with 
charged tracks. Resolved tt^ s are reconstructed by combining pairs of photon candidates and 
requiring that the invariant mass of the resultant candidate is between lOOMeV/c^ and 170MeV/c^. 
Photon candidates used in vr*^ reconstruction are required to have a minimum energy of 50MeV. 
In addition to resolved tt*^ s, unresolved vr" s are identified as single clusters of energy deposited in 
the EMC without distinct local maxima and satisfying certain topological criteria . For the final 
event selection, we tighten the requirement on the vr'' mass to between 120MeV/c2 and ISOMeV/c^. 

Neutral kaons are reconstructed through the decay — > 7r"^7r~. candidates are formed 
from combinations of two oppositely charged tracks satisfying basic track criteria similar to those 
mentioned above, except with a looser constraint on their proximity to the beam spot. We require 
that the candidate's flight length exceed 2 mm, the angle between its flight and momentum 
directions be less than 40 mr, and the mass lie within iblOMeV/c^ of nominal. 

We reconstruct uj mesons using the dominant decay channel, u) — > 7r~7r"'"7r*^, which has a branch- 
ing fraction of 88.8%. Candidates are obtained as combinations of two charged tracks with oppo- 
site signs and one vr*^ candidate. The invariant mass of the uj candidate is required to be within 
0.05 GeV/c^ of the known UJ mass H. Figure H shows the vr vr^vr invariant mass for uj candi- 
dates with a center-of-mass momentum between 1.9 and 3.1GeV/c. The signal is fitted with a 
Breit-Wigner function with the natural width of the u) (8.4MeV/c^) convoluted with a Gaussian 
resolution function, and the combinatorial background is fitted with a second-order polynomial. 
The Gaussian resolution function is found to have a width of 8MeV/c^. 

Candidate p and K* mesons are reconstructed by combining pairs of appropriately charged 
tracks and/or vr*^ candidates, with the assumption of the relevant final-state rest masses, and 
requiring that the combination have an invariant mass sufficiently close to the known resonance 
mass. Requirements for the invariant masses are analysis-dependent, and are described below. 

The r] meson is reconstructed by combining pairs of photon candidates with a minimum energy 
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Figure 1: Mass spectrum of hard to candidates in on-resonance data, fitted with a Breit-Wigner 
convoluted with a Gaussian for the signal and a second-order polynomial for the background. 



of 50 MeV. The reconstructed r] mass for data is shown in Fig. ^. We obtain a width of 19 MeV/ 
for a Gaussian fit, with the background described by a second-order polynomial. The rj' meson is 
reconstructed in the r/' rjir'^TT" decay mode. This is done by combining rj candidates with pairs 
of oppositely charged tracks. A mass constraint is used to improve the rj mass resolution. Fig. |^ 
shows the invariant mass distribution for all r/vrvr candidates with momentum above 2 GeV/ c. From 
a fit to the sum of a Gaussian and a second-order polynomial, we obtain a width of 4.7 MeV/ for 
the rj' peak. 

For B decays to final states formed from a pseudoscalar and a vector meson, the vector meson 
is polarized. We make use of the angular distribution of the vector meson decay products to 
distinguish between signal and background. We compute in the rest frame of the vector meson the 
cosine of the angle Oh between the direction of one of the vector meson daughters (the normal to the 
decay plane in the case of the to) and the direction of the B meson in the same frame. We expect 
signal events to be distributed according to cos^ 9h, while background should be approximately 
flat. 

To veto electrons, we require that charged hadrons satisfy E/p < 0.9, where E is the energy 
measured in the calorimeter associated with a charged track of momentum p. We also require that 
kaons be positively identified as such with a combination of the information from the DIRC and 
the drift chamber dE/dx. For pions, we simply require that they not be identified as kaons. 

Reconstruction of B candidates is done by forming all combinations of the appropriate final-state 
candidate particles, and requiring them to satisfy kinematic constraints appropriate for B mesons. 
We use two kinematic variables for this: mss = \J {^s + Pq • pBy/E^ — p^, where the subscripts 
and B refer to the e~^e~ system and the B candidate, respectively; and AE = E'^ — y/s/2, where 
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Figure 2: Distribution of (a) 77 {a = 19 MeV/c^) mass, and (b) tjtt'^tt {a = 4.7 MeV/c^) mass in 
on-resonance data. 

-E^ is the B candidate energy in the center-of-mass frame. For signal events, the former has a value 
close to the B meson mass Q and the latter should be close to zero. In our analyses, the appropriate 
final state masses are assigned to the final state particles in calculating AE. 

We have measured the expected distributions of the signal with respect to tties and AE, using 
two calibration modes described in Section ^. Based on these, we define a "signal region" in the 
niES-^E plane by |mES- ("ies) | < eMeV/c^ and \ AE-{AE) \ < 70MeV, where (rriEs) and (AE) 
indicate the mean values obtained in the calibration modes. This signal region is common to all 



modes, except B~^ — > ujh'^ (see Section 6.2). 



4 Calibration with charmed B decays 

The decay mode B+ D^tt+,D° K+n' results in a final state containing three charged tracks, 
while the same B+ decay with K+ vr TT yields in addition a vr . The D mass is small 

compared to that of the B, so that this decay has kinematics similar to charmless B decays to two 



^ A nearly equivalent constrained mass that we use for some of the modes is ttiec = \/ Eg^ — Pg^, obtained from 
a fit constrained by Eg = i-\/s- 



11 



light mesons such as B ^ pn, B — > K*'k or i? — > (jtt. The product branching fractions are known 
to about 10% precision and are roughly an order of magnitude larger than those expected of many 
charmless modes. A study of these decays can therefore be used to calibrate the variables used in 
the charmless analyses. The most significant difference is that the has a measurable lifetime, 
resulting in decay tracks which do not share a common vertex with the pion from the B decay. 

The candidate selection proceeds along similar lines to those described earlier for the charmless 
decay modes, except that the invariant mass of the daughters, assumed to be ttK{j:^), must be 
consistent with the mass. The background is estimated by counting events in a "sideband" of 
the m-E,s~^E distribution, and also by looking at K-ktt{i:^) track combinations which satisfy the 
selection criteria but which have an incorrect charge combination (the kaon charge is opposite to 
that of the B meson). The extrapolation from sideband to signal region is done in the same way 
as described for the charmless analyses. 

The efficiency of the selection for B+ D^-k+,W K+Ti~[K^] is (23 ± 2)% [(12 ± 1)%], 
evaluated as described in Section |6|. When combined with the total number of B decays the product 
branching fraction for B'^ D^tt^,D^ K~^7r~ -[vrO] is measured to be (2.21 ± 0.11) x 10"^ 
[(6.79 ± 0.32) X 10"^], where the errors are statistical only. These are both consistent with the 
corresponding world-average values [^ . 

The selected signals in these modes have been used to estimate the rriEs and AE resolution 
functions of our target modes in data. We find the m^s resolution to be 2.7 it 0.1 MeV/c^ in the 
mode with only three charged tracks and 3.0ib0.1 MeV/c^ in the mode with a vr*^. The corresponding 
AE resolutions are 26 it 2 MeV in the mode with only three charged tracks and 33 it 2 MeV in the 
mode with a vr*^. The measured values were used to inform the choice of signal region, noted at the 
end of Section |3| and for calculation of the selection efficiency and its uncertainty. 

5 Background characterization and suppression 

Charmless hadronic modes suffer very large amounts of background from random combinations 
of tracks, mostly from light quark and charm continuum production. Such backgrounds may be 
reduced by selection requirements on the event topologies computed in the T(45) rest frame. We 
use the angle 9t between the thrust axis of the B meson decay and the thrust axis of the rest of 
the event. For continuum-related backgrounds, these two directions tend to be aligned because the 
reconstructed B candidate daughters generally lie in the same jets as those in the rest of the event. 
By contrast, in B events, the decay products from one B meson are independent of those in the 
other, making the distribution of this angle isotropic. In consequence, requiring that this opening 
angle be significant provides a strong suppression of continuum backgrounds. 

Other event shape variables also help to separate signal and background. We combine several 
variables in a Fisher discriminant J- [^. The variables contained in are: 

• The summed energy in nine annular cones of angular extent 175 mrad, coaxial with the thrust 
axis of the B candidate. 

• The absolute value of the cosine of the angle between the B direction and the beam axis. 

• The absolute value of the cosine of the angle between the thrust axis of the B candidate and 
the beam axis. 

The Fisher discriminant is a linear combination of these 11 variables. The coefficients for each 
variable are chosen to maximize the separation between training samples of signal and background 
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events. These coefficients were determined for the i?"*" — > ujh'^ analysis and the same coefficients 
are used for the other oj and r]' analyses. 

Despite the power of such topological variables to reduce the combinatorial backgrounds, most 
of the modes we have searched for continue to suffer significant background levels. Even after 
stringent selection criteria have been applied, it is necessary to do a background subtraction to 
isolate a signal or set an upper limit. In order to do this, the background in the signal region 
must be estimated. This is done by noting that the amount of background in the signal region 
should be related to the amount in a sideband region, located near the signal region in the mES- 
iS.E plane. We define A, to be the ratio of the number of candidates in the signal region to 
the number in the sideband region. Two different methods have been used to estimate A. In 
the first, the shape of the distribution of the background as a function of tties is measured from 
on-resonance data. This is done with events slightly displaced from the signal region in the IS.E 
variable (0.1 < \IS.E \ < 0.3). In the second method, off-resonance data are used, simply counting 
the numbers of candidates in the signal and sideband regions to provide the ratio. The first 
method is used in the three-body measurements, and the second method in the quasi-two-body 
measurements. In each case, the alternative method has been used to evaluate systematic errors. 
Where insufficient data are available for the off-resonance studies, selection criteria in uncorrelated 
variables have been loosened, to provide additional statistics in the determination of A. We find 
that the value of A is quite independent of the decay mode provided the selection criteria are 
common. For some modes with limited statistics we substitute a more precise value of A from a 
mode in which it is better measured. Simulated light-quark (u, d, s, c) continuum events are used 
to check the value of A and provide a statistically-independent sample for optimization of event 
selections in the final stage of the analysis. 



6 Analysis 

The branching fractions are calculated according to 

S=^^ (1) 

where A'^i is the number of candidates in the signal region for on-resonance data; is the number 
of candidates in on-resonance data observed in the sideband region, so that is the estimated 
number of background candidates in the signal region; is the number of BB pairs produced 

and e is the signal efficiency. Except in particular cases, explained below, the sideband is specified 
by 5.20 < mEs < 5.27GeV/c^, \AE - (AE) \ < 0.2 GeV, where (AE) was defined in Section |. 

The numbers and distributions of candidates within the signal region remain unknown to us 
until all aspects of the analysis are finalized. The final selection criteria are chosen to maximize the 
sensitivity of the signal, defined as the expected signal yield divided by its statistical uncertainty. 
Once chosen, neither the background description nor the cuts are changed. This procedure has 
been followed independently for each channel. 

For the signal efficiency in Eq. (|l]) , we use simulated signal events and the same selection criteria 
as used for the data. The efficiencies due to tracking, particle identification and the AE and mEs 
selection criteria are determined from independent control samples from the data. 
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6.1 Modes with three-body final states 

The selection criteria which are varied during optimization are the thrust- angle requirement, the 
helicity- angle requirement and the mass requirement on resonance candidates. In addition to these 
and the selection criteria described in Section ^, a quality requirement of < 20 is applied to the 
three (two) charged track vertices in all charged (neutral) three-body analyses. Finally, a veto on 
all pairs of tracks which are consistent with the — > K^7r~ decay hypothesis (independent of 
particle identification) is applied in all three-body analyses. This has no significant effect on our 
results in any mode except — > p^K^ . The results for these analyses are summarized in Table |l]. 
The branching fraction central values are calculated from Eq. |l[ The significance values given are 
the probabilities of a fluctuation of the background to account for the observed yield, in equivalent 
Gaussian standard deviations. 

The background extrapolation is validated, albeit with low statistics, by the observed ratios 
of events in signal and sideband regions in off-resonance data. In each case, we found that the 
number of events falling in the signal region was consistent within one standard deviation with the 
expectation from the sideband regions. 

The p^K^ mode has the additional complication of cross-talk background from the 

— > r]'K~^ mode with 7]' p^^. The 7 can be quite soft and in such cases the p*^ and 
may imitate the B~^ p^K~^ signal. We correct for this by processing simulated events from this 
background mode using the full set of optimized selection criteria for the B~^ — > p^K~^ mode. We 
obtain an efficiency for this mode of 2.4%, leading to a reduction of the observed branching fraction 
by 2.0 X 10^^. The existence of this S-related background was foreseen, and for this reason the 
sideband used to estimate the background is somewhat smaller than that used for the other modes, 
to ensure that there was no significant leakage of this background into the light-quark background 
estimate. Cross-talk among the other modes considered is found from Monte Carlo simulation to 
be negligible at the current level of precision. 

The distributions of mss a-nd AE for the resonant three-body modes are shown in Fig. ^. 

We note that the analyses of the two modes B~^ — > K~^tt~tt~^ and B~^ 7r'''7r~7r''' differ in some 
respects from those with a resonance. We remove all two-body combinations with invariant masses 
less than 2 GeV/c^, and in addition veto combinations of tracks consistent with the decay mode 
B^ — > J/ipK~^. The distributions of rriES and for the non-resonant modes are shown in Fig. |^. 

6.2 Modes with uj and r]' 

We summarize the results of the measurements in Table |2[ For all decay modes, we require 
I cos 0^1 < 0.9, and then find the optimum selection requirement in JF. The latter is listed for 
each mode in the table. Also listed are the selection criteria on resonance masses and vr*^ or rj 
masses. The selection criteria are centered on the known masses @. We also require that | cos 9h \ 
be greater than 0.55 (0.4) for the decay modes B~^ — > u!h~^{B^ — > loK^). 

Table |^ gives the yields in the signal and sideband regions, both for on- and off-resonance 
data. The signal region is defined as |mES - 5.279| < 0.006 GeV/c^ and |A^| < 0.070 GeV, 
except for 5+ ujh^ where we use -0.113 < A.E < 0.070 GeV, allowing for the 43 MeV shift of 
B~^ — > ijK~^ when the energy of the is computed with a pion mass. The sideband is defined 
as 5.20 < niES < 5.27 GeV/c^ and \AE\ < 0.2 GeV. The branching fraction central values and 
significance values are calculated in the same way as for Table |. 

In Fig. ^ we illustrate the reconstruction of the B~^ r]'K~^ states with plots of the kinematic 
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Table 1: Results for the three-body final state analyses. Energies and masses are in units of MeV. See the text for an explanation 
of the quantities given. 



Quantity B+ K*°7r+ B+ p^K+ B+ K+tt-tt+ B+ p^ir^ B+ 7r+7r-7r+ B^ p^n^ 



cos 6t 


0.70 


0.60 


0.70 


0.55 


0.70 


0.50 


K*^ / p mass 


±100 


±200 


N/A 


±200 


N/A 


±150 


K*^/p helicity angle 


0.40 


0.40 


N/A 


0.40 


N/A 


0.30 


Signal events 














On-vp*^ data 


22 


25 


33 


64 


32 


77 


Off-res data 


1 


1 


3 


4 


5 


6 


Sideband events 














On-res data 


319 


256 


454 


1057 


722 


1129 


Off-res data 


50 


45 


60 


166 


92 


183 


A 


0.037 


0.056 


0.037 


0.037 


0.037 


0.037 


Est. BG 


11.8 ±0.9 


14.3 ± 1.2 


16.7± 1.1 


39.1 ± 1.9 


26.6 ± 1.4 


41.5 ±4.3 


Signal 


10.2 ±4.8 


10.7 ±5.1 


16.3 ± 5.8 


24.9 ±8.2 


5.4 ± 5.7 


35.5 ±9.8 


e X secondary Bi (%) 


9.7 ±0.9 


9.8 ± 1.0 


6.0 ±0.6 


11.9 ± 1.2 


8.2 ±0.8 


8.3t\i 


Stat. sign, (a) 


2.4 


2.2 


3.2 


3.3 


0.7 


4.5 


cross-talk 














correction (xlO~^) 




-2.0 












13±6±1 


10±6±2 


31 ±11 ±3 


24 ± 8 ± 3 


7.5tli ± 0.8 


49 ± 13t^ 


90% CL hmit 


< 28 


< 29 


< 54 


< 39 


< 22 
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Figure 3: Distributions in tties and AE after all other selection criteria have been applied for 
(a, b) B+ K*^7r+ , (c, d) B+ p^K+ , (e, f) B+ pO^+ , (g, h) 5° ^ p^7r± . 
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Table 2: Results for the quasi-two-body decay analyses. Energies and masses are in units of MeV. 
See the text for an explanation of the quantities given. 



Quantity 


B+ ujh+ 




B+ -> rj'K+ 




T 


-0.4 


-0.2 


0.3 


0.6 


uj/-q' mass 


±20 


±20 


±10 


±10 


77° /?7 mass 


±15 


±15 


±45 


±10 


Signal events 










On-res data 


13 





14 


2 


Off-res data 


1 











Sideband events 










On-res data 


128 


19 


45 


14 


Off-res data 


16 


6 


8 


1 


A 


0.055 


0.042 


0.042 


0.042 


Est. BG 


7.1 ±0.9 


0.8 ±0.2 


1.9 ±0.3 


0.6 ±0.2 


Signal 


5.9 ±3.6 


-0.8 ±0.0 


12.1 ±3.7 


1.4± 1.4 


Eff. e (%) 


8.5 


7.6 


17.1 


10.0 


Secondary B (%) 


88.8 


30.5 


17.2 


6.0 


e X secondary B {%) 


7.5 ± 1.4 


2.3 ±0.4 


2.9 ±0.4 


0.60 ±0.09 


Stat. sign, (cr) 


1.7 


0.0 


5.3 


1.1 


6(xlO-6) 


8.9 ±5.4 ±2.2 


0.0 


62 ± 18 ± 8 


27 ±27 ±5 


90% CL limit 


< 24 


< 14 




< 112 
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Figure 4: rriES and AE distributions for (a, b) — > vr^vr vr"^ and (c, d) — > K^tt tt^ non- 
resonant analyses after the final selection. 
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Figure 5: Kinematics of B'^ r]'K^ , rj rjTr^ir : (a) AE vs. m-ES and (b) 77 vs. jyvr'^vr invariant 
mass. 



variables associated with the meson and daughter resonance. For the other decay modes we 
report our final results as 90% confidence upper limits. 

The distributions for all decay modes, projected on to the rriES and AE axes, are presented in 
Fig. |. 

7 Systematic errors 

The primary uncertainty in the signal yield is the counting statistics of the accepted events in the 
signal region. The systematic uncertainty on the background estimate receives contributions from 
the number of events observed in the on-resonance sideband region and from the ratio A. The 
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uncertainty on the latter depends on the statistics of the off-resonance samples and on the method 
used to determine it. We estimate an overall uncertainty in A of 4% for the charged three-body 
modes and 10% for the other modes. 

The branching fraction itself has a multiplicative uncertainty from our knowledge of the effi- 
ciency, e. We measure e using signal Monte Carlo simulation, and a significant uncertainty arises 
from limited signal mode statistics, the fractional error due to the size of the Monte Carlo samples 
used to estimate e varies between 4% and 7%, depending on the mode. The accuracy of the sim- 
ulation is subject to systematic uncertainties in the efficiency of tracking and calorimetric shower 
detection, particle identification criteria, and the resolutions on quantities used to define selection 
criteria. At present the dominant uncertainties are in the tracking, tt^ /r] reconstruction, and Fisher 
Discriminant selection efficiency, with some variation on their relative importance depending on 
the decay mode. 

Independent studies determine an absolute track-finding efficiency uncertainty of 2.5% per track, 
added coherently for all tracks required in the analyses. For vr'^ or r] reconstruction, we assign an 
uncertainty of 5%, and add an additional uncertainty of 1-5% in quadrature to account for the 
effect of the final mass requirement. For the Fisher Discriminant selection, the uncertainty on 
the efficiency can vary considerably, depending on the tightness of the selection requirement. The 
uncertainty ranges from only a few percent for decay modes with low background, to as much as 
15% for ones which need a tight selection requirement. The number of produced BB events has 
been estimated in separate studies with an estimated uncertainty of 3.6%. Uncertainties in 
the particle identification and uies ^-nd AE selection criteria contribute systematic uncertainties 
of 2-4%, depending on mode. The overall systematic uncertainty is the quadrature sum of the 
contributions from all sources. 

For modes in which we find an upper limit for the branching fraction, we account for the 
uncertainty of the efficiency by using in Eq. |^ the measured efficiency reduced by one standard 
deviation in the systematic error. 

8 Summary 

We have presented a number of preliminary measurements of charmless hadronic B decays, summa- 
rized in Table |^. We observe significant signals for B~^ rj'K^ and B^ P^tt^, with branching 
fractions B{B+ r]'K+) = (62 ± 18 ± 8) x 10^^ and B{B'^ p^n"^) = (49 ± 13^1) x 10"^. Upper 
limits are given for the remaining channels we have studied. In each case, our determinations are 
consistent with earlier published measurements |Q, H, |6|. 
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Table 3: Summary of branching fraction measurements. Inequality denotes 90% CL upper limit, 
including systematic uncertainties. 
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